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SUMMARY 

A general calculation is given for the earth’s albedo input to 
a spherical satellite, with the assumption that the earth can be 
considered a diffusely reflecting sphere. The results are presented 
in general form so that appropriate values for the solar constant 
and albedo of the earth can be used as more accurate values become 
available. The results are also presented graphically; the incident 
power is determined on the assumption that the mean solar constant 
is 1.353 x 10 6 erg cm -2 sec -1 and the albedo of the earth is 0.34. 
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EARTH REFLECTED SOLAR RADIATION 
INPUT TO SPHERICAL SATELLITES 

by 

F. G. Cunningham 

Goddard Space Flight Center 


INTRODUCTION 


Several recent papers have considered the thermal behavior of satellites and space 
vehicles in near-earth orbits. Goldman and Singer, considering the temperature environ- 
ment of a spherical satellite, made the approximation that a satellite receives reflected 
solar radiation from only that portion of the earth directly beneath it (Reference 1). This 
does not take into account the diffusiveness of reflection from the earth. Wood and Carter 
calculated the reflected solar radiation input to a sphere located on the earth- sun line, 
but not for any other orientations (Reference 2). Hibbs’ paper (Reference 3), which includes 
the basic integral expression for reflected solar radiation, gives a result which is errone- 
ous. This error will be reviewed presently. 


CALCULATIONS 

The geometry of the problem is shown in 
Figure 1. The associated definitions are: 

r = The distance of the satellite from the 
center of the earth, in mean earth radii. 

p ~ The distance from the satellite to the 
element of area ds on the earth’s sur- 
face, in mean earth radii. 

9 = The co- latitudinal angle defining the 
position of ds with respect to r. 

<t> = The azimuthal angle of integration de- 
fining ds with respect to the plane 
formed by r and the solar vector. 

€ = The angle between p and the normal to 
the surface element ds. 



Figure 1 - Geometry of the problem 
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k = The angle between the solar vector and the normal to the surface element ds. 

8 = The angle between the solar vector and r. 

S = The solar constant. 
a = The mean albedo of the earth. 

The general expression for the reflected solar radiation incident upon a satellite, p se , can 
be written as 


^ Sa j cos k cos £ ds 

*SE “ IT J p2 


( 1 ) 


where the cross sectional area of the satellite is assumed to be unity. From Figure 1, 
the quantities in Equation 1 are: 

p 2 = r 2 + 1 - 2r cos 8 ; 


ds = sin 8 dd dcp 


( 2 ) 


cos £ 


r cos <9-1 

~ j 

(r 2 + 1 - 2r cos 8)^ 


and 


cos k = cos 8 cos <9g + sin 6 sin <9 g cos <p 


(3) 


Equation 3 follows from the relationships of the angles and sides of a spherical triangle. 
Substituting Equations 2 and 3 into Equation 1 gives the complete expression for the re- 
flected solar radiation, 


SE 




(r cos 8-1 ) (cos 6 cos 8 S + sin 8 sin 6 s cos <p) sin 0 d 6 d<p 

3 

(r 2 + 1 - 2r cos 8) 2 


(4) 


Equation 4 is to be integrated over the area of all the terrestrial surface elements 
which both "see" the satellite and are illuminated by the sun. The limits of integration for 
8 and cp in Equation 4 depend on the relations between the two relevant regions: Region 1, 
the hemisphere of the earth that is illuminated by the sun; and Region 2, the spherical cap, 
consisting of all the elements of the earth which see the satellite. All elements of Region 
1 satisfy the condition k <tt/2 . The great circle which constitutes the boundary between 
the sunlit and unlit hemispheres, i.e., the terminator, is the locus of points for which 
k = 77/2 • All elements of Region 2 satisfy the condition 0 < 6 < 8 m , where 9 m = cos -1 l/r . 

There are three possible relations between Regions 1 and 2 which define the areas of 
integration: 

a. Region 2 lies completely within Region 1; 

b. Region 2 lies partly within Region 1; or 

c. Region 2 lies completely outside of Region 1. 
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For condition a, o < # s < V2- 0 m , and takes on all values from o to 2 tt ; but because 
of symmetry with respect to the plane determined by the earth- satellite and earth- sun lines, 
the integration can be performed between the limits o and n and the result multiplied by 2 . 
In this case o < 6 < & m . In condition c, where 0 g > 77/2 + 6 m , no reflected sunlight is inci- 
dent on the satellite; and this condition is of no further interest. 


Condition b presents two possibilities: b 1? in which 77/2 - d m < d s < V 2 ; and b 2 , in 
which 77/2 < 0 S < V2 + 0 m . For b x the area of integration is divided into two parts: b n for 
which 0 i 6 < V 2 - # s and 0 < <£ < 77; and b 12 for which 77/2 - 6 S < e < 6 m . In part b 12 of 
subregion b 1? and in subregion b 2 where d s - n/2 < 6 < e m , the upper limit of the <f> inte- 
gration, c£ m , may be determined by solving Equation 3 for cos 4 > m after substituting the 
value ^ = 77/2 : 


cos & 

1 m 


cos 6 cos # s 
sin d sin < 9 g 


( 5 ) 


The fact that the range of <fi is not 0 to 277 (that is, 0 to 77 in the present paper) for all 
values of 6 in the range 0 < 6 < d m was overlooked in the paper by Hibbs (Reference 3 ). The 
result, upon integration, is that the term in Equation 4 containing cos 4 > is eliminated for 
all values of e s . Consider Equation 4 when - 77/2 • If the cos <£ term is zero upon inte- 
gration and the cos $ s term equals zero by definition of the problem, then the incident power 
is identically zero. But, clearly, when e s = 77/2 , half of Region 2 lies within Region 1 and 
the satellite does indeed see illuminated portions of the earth. 


In summary, the regions and subregions of integration are: 

77 0 < 6 < e m 

a. 0 1 ^ 1 r - 0 m , for which ’ 

S 2 m ’ O<0<77, 


and the incident power = P gE 


( 6 ) 


bj- ^ < 6 S <^ , where bj = b u + b J2 , (7) 

bjj^. 0 < 6 < ^ - and 0 <. <f> < 77 , 

b 12‘ \ ~ and 0 < 4 > < cos -1 (-cot 6 cot < 9 S ) , 

and the incident power = P gE ; 


2 ' 




JL ^ a 

2 + ° n 


for which 


e s ~ ~ < e < e m 

0 < <£ < cos“* (-cot 6 cot ( 9 g) 


and the incident power = P gE 


( 8 ) 
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In consideration of the foregoing, and with the inclusion of Equation 3, Equation 4 yields 


SE, 7 7 


■^f f 

J o J o 


(r cos 0-1) (cos 6 cos 0 S + sin 6 sin 0 S cos <£) sin 6 d0 d 4> 

— j- 

(r 2 + 1 - 2r cos d) 2 


( 9 ) 


where 


3 = 2S< 

SE- “ 77 


C 

Jr, Jr 


0 < e s 


(r cos 6-1) (cos 6 cos % + sin 6 sin 0s cos 4>) sin 6 d6 d 4> 


(r 2 + 1 - 2r cos 6 ) 


2 S c 


where 

and 

p - 

r SE„ " 77 


/» 0m -*1411 

0 


(r cos 6-1) (cos 6 cos 0g + sin 6 sin 6 g cos <£) sin 6 d6 dfi 

— 

(r 2 + 1 - 2r cos 0) 2 


( 10 ) 


— - e m < 6c < 


r 0 m J- 

L •’o 


s T ' 


(r cos 0 - 1) (cos 0 cos 0g + sin 0 sin cos <£) sin 6 dd dcfi 


— < 6 S < — + 6„ 

2 S ~ 2 " 


b - * 

where 

Upon integration, Equations 9, 10, and 11 become: 

( 2r+ ^)- ( 2+ -^)( r2 - 1 )"] cos ‘ ?l 


(r 2 + 1 - 2r cos 6 ) 


2§c 


s ; 


( 11 ) 


( 12 ) 


2§a cos 0s 


•SE, 


^1 + -i- j ^-2r 2 + 1 j ^ 2r sin 0 S + sin 2 0 S ( 

■ ~j~ 

(r 2 + 1 - 2r sin 0g) 2 


+ — + 2r 


2$a f 

77 J/. 


(r cos 0-1) cos 0 S sin 0 cos 0 / cos 0 cos 0 S 


— - 0 . 


(r 2 + 1 - 2r cos 0) 


J / cos 0 cos 0g \ 
l sin 0 sin 0s ) 


dd 


2§a 


1 


— - 0 * 
2 


(r cos 0-1) sin 0 (sin 2 0 - cos 2 0s) 2 

- ;t 

(r 2 + 1 — 2r cos 0) 2 


d0 ; 


(13) 
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m 


(r cos (9-1) cos (9g sin 9 cos 6 

— 

(r 2 +1 - 2r cos o ) 2 


sm U sin 


% -i) 


(r cos 0-1) sin Q (sin 2 0 - cos 2 0 S ) 
" T 

(r 2 + 1 - 2r cos &) 2 


SIMPLIFICATION OF THE EQUATIONS 


Equations 12, 13, and 14 may be simplified by defining the following quantities: 


— cos <9 S 


-2r 2 + 1] + (2r - — | sin B c + sin 2 0, 


(r 2 + 1 - 2 


r sin # s ) 2 


■-1 


(r cos 6-1) cos 9 S sin 9 cos 9 , 

— cos 1 


^ 2 — (r 2 + 1 - 2r cos 6 ) 2 


cos 9 cos (9g 
sin 9 sin 9 S 


d 6 , 


*1 




(r cos (9 - 1) sin 6 (sin 2 9 - cos 2 6$) 2 

3 

(r 2 + 1 - 2r cos 9] 2 


d6 , 


where 9 S has the range indicated in Equation 7; 


(%-t) 


(r cos 6-1) cos 6~ sin 6 cos 9 « / cos 6 cos 6c \ 

— ° cos" 1 - L) d 6 , 

3 y sm 9 sin 0g J 

(r 2 + 1 - 2r cos o) 2 


m (r cos 9-1) sin 6 (sin 2 6 - cos 2 ( 9g) 2 d<9 

&S ~~%J (r 2 + 1 - 2r cos 6) 2 
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where 0 S now has the range indicated in Equation 8; and 

i- J_n * 



where the range of d s is given in Equation 6. 

With these definitions, the inputs can be written simply: 


P SE X = SaF 5 

(21) 

= Sa(A + B + C) ; 

(22) 


and 


p se 3 = Sa ( D + £ ) • ( 23 ) 

The error introduced by the assumption of a perfectly spherical earth should be small and 
is neglected. However, the value used for the average albedo (a) can introduce a sizeable 
variation into the result. The hemispherical yearly average of the earth 1 s albedo seems to 
vary between 0.34 and 0.40, depending upon the investigator (Reference 4). 


RESULTS 


Figures 2 and 3 show the results of the preceding calculations, with the value 0.34 for 
the average albedo of the earth and the value 0.1353 watt/cm^ for the solar constant. The 
power input (for unit absorptivity) per unit cross-sectional area of the satellite is plotted 
as a function of altitude above the surface of the earth. The angle 3 S is the parameter 
which generates the family of curves. Figure 2 presents the results in the range from 
200 to 3200 kilometers and Figure 3 covers the range from 200 to 32,000 kilometers. 
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ALTITUDE (102 KILOMETERS) 

Figure 2 - Power input per unit cross-sectional area of a satellite plotted as a function of altitude 

for the range 200—3200 miles 
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Figure 3 - Power input per unit cross-sectional area of a satellite plotted as a function of altitude 

for the range 200—32,000 miles 

REFERENCES 

1. Goldman, D. T., and Singer, S. F., "Studies of a Minimum Orbital Unmanned Satellite 

of the Earth; Part IV, Radiation Equilibrium and Temperature", Physics Dept. Technical 
Report No. 46, University of Maryland, College Park, Maryland, June 1956 

2. Wood, G. P., and Carter, A. F., "Predicted Characteristics of An Inflatable Aluminized- 
Plastic Spherical Earth Satellite With Regard to Temperature, Visibility, Reflection 

of Radar Waves, and Protection from Ultraviolet Radiation", NASA Technical Note 
D-115, October 1959, pp. 12 ff 


D-iuyy 


D-1099 


9 


3. Hibbs, A. R., ’’The Temperature of an Orbiting Missile”, Jet Propulsion Laboratory, 
Pasadena, California, Progress Report No. 20-294, March 28, 1956, p. 3 

4. Johnson, J. C., ’’Physical Meteorology", New York: Wiley, 1954, p. 171 


NASA-Langley, 1961 


NASA TN D-1099 

National Aeronautics and Space Administration. 
EARTH REFLECTED SOLAR RADIATION INPUT TO 
SPHERICAL SATELLITES. F. G. Cunningham. 
October 1961. 9p. OTS price, $0.50. 

(NASA TECHNICAL NOTE D-1099) 

A general calculation is given for the earth’s albedo 
input to a spherical satellite, with the assumption that 
the earth can be considered a diffusely reflecting 
sphere. The results are presented in general form 
so that appropriate values for the solar constant and 
albedo of the earth can be used as more accurate 
values become available. The results are also 
presented graphically; the incident power is deter- 
mined on the assumption that the mean solar constant 
is 1.353 x 10® erg cm“2 sec"l and the albedo of the 
earth is 0.34. 


I. Cunningham, F. G. 

II. NASA TN D-1099 

(Initial NASA distribution: 
7, Astrophysics; 

21, Geophysics and 
geodesy; 26, Materials, 
other; 35, Power sources, 
supplementary; 

47, Satellites; 

52, Structures.) 


NASA TN D-1099 

National Aeronautics and Space Administration. 
EARTH REFLECTED SOLAR RADIATION INPUT TO 
SPHERICAL SATELLITES. F. G. Cunningham. 
October 1961. 9p. OTS price, $0.50. 

(NASA TECHNICAL NOTE D-1099) 

A general calculation is given for the earth’s albedo 
input to a spherical satellite, with the assumption that 
the earth can be considered a diffusely reflecting 
sphere. The results are presented in general form 
so that appropriate values for the solar constant and 
albedo of the earth can be used as more accurate 
values become available. The results are also 
presented graphically; the incident power is deter- 
mined on the assumption that the mean solar constant 
is 1.353 x 10® erg cm~2 sec - l and the albedo of the 
earth is 0.34. 


I. Cunningham, F. G. 

H. NASA TN D-1099 

(Initial NASA distribution: 
7, Astrophysics; 

21, Geophysics and 
geodesy; 26, Materials, 
other; 35, Power sources, 
supplementary; 

47, Satellites; 

52, Structures.) 


Copies obtainable from NASA , Washington 


NASA 


Copies obtainable from NASA , Washington 


NASA 


NASA TN D-1099 

National Aeronautics and Space Administration. 
EARTH REFLECTED SOLAR RADIATION INPUT TO 
SPHERICAL SATELLITES. F. G. Cunningham. 
October 1961. 9p. OTS price, $0.50. 

(NASA TECHNICAL NOTE D-1099) 

A general calculation is given for the earth’s albedo 
input to a spherical satellite, with the assumption that 
the earth can be considered a diffusely reflecting 
sphere. The results are presented in general form 
so that appropriate values for the solar constant and 
albedo of the earth can be used as more accurate 
values become available. The results are also 
presented graphically; the incident power is deter- 
mined on the assumption that the mean solar constant 
is 1.353 x 10® erg cm"2 S ec“l and the albedo of the 
earth is 0.34. 


I. Cunningham, F. G. 

II. NASA TN D-1099 

(Initial NASA distribution: 
7, Astrophysics; 

21, Geophysics and 
geodesy; 26, Materials, 
other; 35, Power sources, 
supplementary; 

47, Satellites; 

52, Structures.) 


NASA TN D-1099 

National Aeronautics and Space Administration. 
EARTH REFLECTED SOLAR RADIATION INPUT TO 
SPHERICAL SATELLITES. F. G. Cunningham. 
October 1961. 9p. OTS price, $0.50. 

(NASA TECHNICAL NOTE D-1099) 

A general calculation is given for the earth’s albedo 
input to a spherical satellite, with the assumption that 
the earth can be considered a diffusely reflecting 
sphere. The results are presented in general form 
so that appropriate values for the solar constant and 
albedo of the earth can be used as more accurate 
values become available. The results are also 
presented graphically; the incident power is deter- 
mined on the assumption that the mean solar constant 
is 1.353 x 10® erg cm"2 sec‘1 and the albedo of the 
earth is 0.34. 


I. Cunningham, F. G. 

II. NASA TN D-1099 

(Initial NASA distribution: 
7, Astrophysics; 

21, Geophysics and 
geodesy; 26, Materials, 
other; 35, Power sources, 
supplementary; 

47, Satellites; 

52, Structures. ) 


Copies obtainable from NASA, Washington 


NASA 


Copies obtainable from NASA , Washington 


NASA 






